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Abstract 

We describe some new exact solutions for two- and four-level systems. 
In all the cases, external fields have a restricted behavior in time. First, we 
consider a method to construct new solutions for one-spin equation and 
give some explicit examples, one of them is in a external magnetic field 
that acts during a finite time interval. Then we show how these solutions 
can be used to solve the two-spin equation problem. A solution for two 
interacting spins is analized in the case when the field difference between 
the external fields in each spin vary adiabatically, vanishing on the time 
infinity. The latter system can be identified with a quantum gate realized 
by two coupled quantum dots. The probability of the Swap operation for 
such a gate can be explicitly expressed in terms of special functions. Using 
the obtained expressions, we construct plots for the Swap operation for 
some parameters of the external magnetic field and interaction function. 



1 Introduction 

Finite-level systems have always played an important role in quantum physics. 
In particular, two-level systems possess a wide range of applications, for exam- 
ple, in the semi-classical theory of laser beams [I], optical resonance [2], and 
nuclear induction experiments [3], and so on. The best known physical system 
that could be identified with a two-level system is a fixed spin-one-half object 
interacting with a magnetic field. The two-level system have been studied by 
many authors using different methods, see for example [H[5]. Likewise the four- 
level systems can be used to describe two interacting one-half spins, e.g., the 
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valence electrons in two coupled semiconductor quantum dots [§J. The most 
detailed theoretical study of the quantum mechanical equations for two and 
four-level systems, and their exact solutions, are presented in [71 [8]. Recently, 
two- and four-level systems have attract even more attention, due to their re- 
lationship to the problem of quantum computation [9] . In this problem, the 
computation is performed by the manipulation of the so-called one- and two- 
qubit gates [12]. The one-qubit gate can be identified with a two-level system 
and two-qubit gates can be identified with a four-level system. For these reasons, 
two- and four-level systems are crucial elements of possible quantum computers. 
For physical applications, it is very important to have explicit exact solutions 
of two- and four-level system equations. In [13] exact solutions of a four-level 
system are used to describe the theoretical construction of a universal quantum 
XOR gate using two-coupled quantum dots. This work shows how the exact 
solutions can be used to establish all the necessary conditions on the external 
fields needed for the implementation of the gate. 

In the present work, we describe some new exact solutions for two- and 
four- level systems that were not represented in our previous works [7J[8]. These 
solutions are found for external fields that have a restricted behavior in time, for 
example, the first solution for two-level system in external field that acts along 
a finite time interval. In Section 2 we describe a general method to construct 
exact solutions for two-level systems with external fields restricted in time, and 
also in this section, we use this method to obtain two explicit external fields 
and the respective exact solutions. In Section 3 we show how these results 
can be applied to construct new exact solutions for the four-level system. This 
system is identified with two interacting spins and we made a detailed study of 
the important case when the difference between the external fields in each spin 
vary adiabatically (vanishes with time). This system can be identified with a 
quantum gate realized by two coupled quantum dots. The probability of the 
Swap operation for such a gate can be explicitly expressed in terms of special 
functions. Using the obtained expressions, we construct plots for the Swap 
operation for some parameters of the external magnetic field and interaction 
functions. 



2 Two-level systems 
2.1 General 

We recall to the reader that two-level systems are described by the so-called 
one-spin equation 

dv 

i— = (a¥)v, «r=(ffi,<T2,o- 3 ), F = (F lt F 2 , F 3 ) , (1) 

where v = v (t) is a two-component spinor, (k = 1, 2, 3) are Pauli matrices, 
Fk = Fk (t) are components of external field strength, see [7j. The general 
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solution of the spin equation reads 

v(t)=R(t)v°, (2) 
where 2x2 matrix R (t) obeys the same spin equation 

.dR(t) 



dt 



(crF)R(t), (3) 



and v° is an arbitrary constant spinor. If R (t = to) = I, where / is 2 x 2 unity 
matrix, then R(t) = ii(t), where u (t) is the evolution operator of the spin 
equation. In the general case, the evolution operator is constructed by the help 
of the matrix R (i) as follows 

u(t) = R (t) R^ 1 , Rq = R(t = t ), u(t = t ) = I . (4) 

The matrix R(t) can always be represented in the form 

R(t) =Ip -i(<rp), p= (Pl,P2,P3), Ps =Ps(t), 5 = 0,1,2,3 , (5) 

The functions p s (t) obey the set of equations 

Po + (pF) = 0, p + [p x F] - p F = 0, (6) 

which follows from ((3j) • Equations ([6]) imply that A = det R (i) = Pq + p 2 is an 
integral of motion. 

Let us suppose that spin equation (J3j) is self-adjoint, which means that the 
external field F is real. In such a case we can chose the functions p s (t) to be 
real. Without loss of generality, in such a case, we can set A = 1, which means 
that R(t) is nonsingular. Under the condition A = 1, the functions p s (£) can 
be expressed via three real parameters a = a (t) , 6 = 6 (t) , and ip — ip(t) as 
follows 



ip — a f 
Po = cos — — cos - , pi 



. ip + a . 8 

sm sm — , 

2 2 ' 



w + a 6 ip — a 9 

P2 = cos — — sin - , p 3 = sin — — cos - . 

With these functions, the evolution operator © assumes the for 



R = R 3 (-<p) R 2 (8) i? 3 (a) , R l [(3 (t)] = exp 



0{t) 



(7) 



(8) 



lr The operators Ri (0) in JSJ are rotation of an angle j3 along the i-axis and the functions 
tp,6 ,a can be identify with the Euler angles. In a similar manner, the functions Pi in JJJ can 
be identify with the Euler parameters. 
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Substituting this expression in J3]) we see that the above R (t) is the evolution 
operator for the spin equation in the following external field 

r, 6 ■ « • A 

rx = X sm f + ~2 sm " cos V > 
d 

i<2 = -jT sm " Sin <y3 — — COS , 

F 3 = |-|cos0. (9) 

For any continuous time-dependent functions <p, 9, a with continuous derivative. 

Of course, the above expressions can not be used to solve the general problem 
for a specific filed, due to the difficulty to solve the integral equations. However, 
these expressions can be used to find exact solutions for external fields with 
some particular characteristic. For example, we can construct exactly solutions 
for periodic external fields by setting 

8 = uit , a = fit , cos ip = C , 

where w, fi and C are constants. Other important kind of external fields are 
those whose action is restricted in time, once these are the most common fields 
in experiments. We can find a variety of solutions for restrict in time external 
fields, just by constructing functions that assume a constant value outside of 
some interval. In the next section we give some explicit realization of this kind 
of fields. 



2.2 Exact solutions for some restricted in time external 
fields 

We can use the results of the above section to construct exact solution of the spin 
equation for external fields restricted in time . It can be done by constructing 
continuous time-dependent functions ip, 9, a, with continuous derivatives, which 
assumes a fixed value outside of a certain interval up = 9 = a — for \t\ > Tj. 

1. Let us chose the functions: 

9(t) = 9 , (p (t) = ip (h a (t) =a , t< -T, 

*(*) = — 2— sm 2T + ^— • I*' <T ' 
V(*) = — 2— sm— + — 2— , \t\<T, 

ai - a . Tit ai + a 
*(*) = — 2— sm— + , \t\<T 

9 (t) = 9x , <p (t) = Vl ,a (t) = a u t > T , (10) 

where cto, a\, 6o, 9i, (po, and ipi are arbitrary constants. With this, the 
external field F |(9|) will be zero at \t\ > T, where T is a constant, and for 
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\t\ < T we have 



h) cos — simp + — (ao - aijcos — sm6»cos^, 

Fa ^ = 8T ^° ~ a i) cos ^ sin6,sin ^- (0Q -0i)cos — cos</>, 
^3 (*) = (Vi - ¥>o) cos — - — (a - ai) cos — cos (9. 

The external field under consideration is not zero only on a finite interval 
\t\ < T and is continuous for all t. The exact solution of the equation (J3j) 
for such a field is constructed by substituting the function (fit)]) in ((7]). 

2. Let the functions 6 — 9(t), tp — <p(t) and a (t) have the form 



a(t) 



Oat 




7^ exp 

J-i 


-(*: 


a t 


"[-(£ 


7^ exp 





3 

n , WW = exp 



(11) 



where ai, 6*o, 0i, fo, fi, ao, and Xfc, (fc = 1,2,3) are arbitrary constants. 
With this, the external field F ((9j) becomes 



*i (*) = - 



F2(t) 



1 - 2 



Ti 

T3 



F 3 (*) = - 
T 2 



T 3 







exp 





1-2 l^r 



1-2 ^ 



exp 



exp 



exp 



exp 



exp 



r 3 
* 

t 

T3 



sin if 



sin cos if 



cos 93 



sm y sin f 



cos 9 



(12) 



This external field vanishes at t — > ±00. The exact solution of the equation 
((3]) for such a field is constructed by substituting the function l[TTj) in ((7]) . 

3. By combining the ip, 9, a functions of the two above examples it is possible 
to construct 6 fields more in the form J9]) that are restricted in time, with 
the exact solutions given by |(5|). 
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3 Four-level systems 



3.1 General 



We write the Schrodinger equation for a four-level system in the following form 
( h = 1), see [8|: 



i^=H(G,F,J)9, 



H = (p-G) 



(S-F) + ^(S.p). 



(13) 



Here ^ is a four-component column; in the general case the interaction function 
J, as well as, the external fields (the vectors G and F) are time-dependent; and 
the 4x4 matrices p and X have the forms 



I ®cr , p = <r(g>/,(X-p)=er(g><r 



3 

£< 

i=l 



<7, . 



where <r = (tri ,0-2,(73) are the Pauli matrices, and / is the 2x2 identity matrices. 
The Hamiltonian matrix reads 



H 



( F 3 + G 



V 



3 "I" ^3 "I" ^ 

Fi + iF 2 
Gi + iG 2 




F 1 - iF 2 

G3 — F3 — 
J 

Gi + iG 2 



Gi - iG 2 
J 

F3 — G3 — 







Gi - iG 2 
Fi - iF 2 
4 - G 3 - F 3 J 



(14) 



Such a model is used to describe two spins subjected to the external magnetic 
fields F and G, and interacting with each other through a spherically symmetric 
Heisenberg interaction whose intensity is given by the interaction function J. 
In particular, this model was used to describe two coupled quantum dots [13] . 
In our work [8] a series of exact solution of equation (fl3"ll for different choices of 
the interaction function and the external fields are found for the first time. 



3.2 Reduction to the two-level system case 

For a special case of two spins subjected to parallel external magnetic fields, 
which we write as 



G = (0, 0, MB9iBi) , F = (0, 0, ^32^2) , B x . 2 = B xa (t) , (15) 

where [Ib is the Bohr magneton and g\ and 172 are g-factors for the corresponding 
spins, one can show that the evolution operator U (t) for the equation (fT3"ll can 
be reduced to an evolution operator u (t) Q for the Schrodinger equation of a 
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two- level system [8]. Such a reduction is given by the equation 



U(t) 

r(t) = 

M = 



exp 
t 



[(S 3 + p 3 )r(t) + S3P3*(i)] M(t) , 



5+ (t) dr , $ (i) 



o 



/ J (i) dr , B + = /i B + g 2 B 2 

Jo 



1 ^ 

Mil Ul2 

M21 U22 

V o o 01/ 



(16) 



where u (t) = \i 
system (see [7]) 



obeys the Schrodinger equation for the following two-level 



flit 

i^ = (*-K)u, fi(0)=I, 



K (t) = (J (t) , 0, S_ (t)) , B_=hb {g 1 B 1 - g 2 B 2 ) 



(17) 



Thus, in the case under consideration, the four-level system problem is reduced 
to solve the two-level system problem lfl7|) with an effective magnetic field K. 

We can now use the expression © to construct exact solution for the four- 
level system l|13p . For an external field K in the form of l|17p we have from 

d ■ a ■ & n 

— sin a sin u> cos tp = (J , 

and a solution in the form (JSj) can be construct for an external field in the form 

V 



K X = J 



2 sin ip ' 



Ko = B_ 



2 tamp 



T) = In (cos ip sin ( 



(18) 



with and a any continuous time-dependent functions with continuous deriva- 
tive. So we can construct restricted in time interactions or fields difference by 
choosing functions 9 and r\ that assume a constant value outside of a desired in- 
terval. Also in this case, explicit examples can be constructed with the functions 
given in the preceding section. 



3.3 An adiabatic variation of the field difference in each 
spin 

Although the expression (fT8"ll allows to construct a variety of external fields 
with some particular characteristic, the solution for a specific field can hardly 
be constructed in this manner. In this section we will analyze a special case for 
a specific external field restricted in time parallel . Consider a four-level system 
in which the field difference varies adiabatically with time, while the interaction 
function is constant. Namely, we chose 

J = a , B- (t) = c/ coshtot , (19) 
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where a, c and lu are real constants. In practical application the pulse applied 
to the system (e.g., two coupled quantum dots) needs to be shorter than the 
decoherence time of the system. But such fast pulse can cause a transition of 
the system to higher energy levels and, consequently, its dynamic can no longer 
be described by the Hamiltonian lfl4|) . The c/ cosh tot dependence is the most 
adequate kind of a variation to avoid this higher energy level transition |14j . 
In addition, it is reasonable to assume that if the only quantity that varies is 
£>-, and B + 3> B-, the interaction function will remain constant [13]. With 
regard to the variation of there are some proposals for the application of 
localized magnetic fields [15] and some techniques that permit the manipulation 
of the g-factor by changing the size of the dots or by the application of external 
electromagnetic fields [16} I17j. 

From the previous Sect., we know that the evolution operator (fl6| of a four- 
level system with the parameters lfl9|) is expressed via an evolution operator of 
a two-level system with effective field 

K (t) = (a, 0, c/ coshwt) . (20) 

The exact solution for the evolution operator with such a field can be constructed 
using our previous results [7]. It has the form 

[ ) ~ |G?| 2 + |G°| 2 I G 2 (z) G 1 (z) J { -Gl G? ) ' {21 > 



where 



Gi (z) = i (2c + lu) z" (1 - z) v F (a, v; 7 ; z) , 

G 2 (z) = 2az^ +1 ' 2 (1 - zf F {a, v + 1; 7 + 1; z) , G° 2 - G x , 2 (-1) , 



\ef - i 



, ip = u>t, a = 7 + v 



M=^,- = ^,7=5+2M. (22) 

F (a, f3, 7, z) is the Gauss hypergeometric function, and complex conjugate quan- 
tities are designated by a bar above. 

Substituting ((22j) into ([TTjl. we obtain 



/ exp[-i(at + T(t))} \ 



uW 

V exp[-i(ot-r(t))] 7 



i?(t) = cx P ( — 



with u (t) given in (f2"Tj) . 

Thus, any transition amplitude for the four-level system can be calculated 
with the help of the evolution operator. Let us, for example, calculate the 
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transition amplitude between the states |t|) and ||t) i which have the form 



|U) = |T>®U> = 



1 

o 



lit) = 11)® IT) = 



o 
i 

V o J 



(23) 



The transition between these states represents, in quantum computation, the 
so called Swap operation and can be experimentally measured [18] , Prom the 
general expression l[T6j) . we see that 



<uiiziiT) = aHi) , uti£iu> = uht) , it) = 



(24) 

Therefore, in the case of the Swap operation between the states l(23|) . we need 
only to calculate matrix elements of the two-level system evolution operator. 
One has also to stress that in this case the Swap operation does not depend on 
the fields' sum B + . 

Using l]2ip. we calculate the probability amplitude for the Swap operation 
with the adiabatic variation lfl9|). 

■ /m-u.v.2 \G 2 (z)G° 1 -G 1 (z)G 2 >\ 2 



(|G°| 2 + |G?f 



In order to use the adiabatic pulse to implement some quantum operations 
(like the Swap or the XOR gate) the duration of the pulse needs to be shorter 
than the dephasing time of the system. For example, in GaAs quantum dots 
this time is about 10 ns [18J, which correspond to u> ~ 1 GHz. In typical exper- 
imental conditions, we have fields of about 5 T, J = 2 x 10~ 3 eV and, to satisfy 
the condition B + 3> B-, we can set the amplitude = 11 mT. So, some 

characteristic values for our system are 



l . 



In Figure 1 we have plots of the probability as a function of time for the 
above values of the parameters. The first maximum occurs at t = 0.5 ns with a 
probability of P = 90%. For larger time, as the cosh -1 approaches to zero, this 
probability varies as A\ sin 2 (at) + A 2 where Aj — Ai(uj,a,c). The amplitude A\ 
decreases as c increases while the shift A 2 increases. The functions A4 change 
significantly with u> only for c > 10a. 
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Figure 1 - Probability of the Swap operation as a function of time for 
J = 2 x 1CT 3 eV, u = 1 GHz and B_ = 11 x 1(T 3 T. 

The dependence of the probability on the parameter u> become noticeable for 
c > 4a. In Figure 2 we plot this dependence for a/c = 2 and a/c = 1/6. The 
parameter u> can be used to significantly attenuate the Swap transition for values 
of c> 10a. 



1.0 - 




Figure 2 - Probability as a function of u> for the values c = 1 (dashed line) and 
c = 12 (solid line) in t = 0.8ns and a = 2. 

A numerical study shows a strong dependence of the maximum values on 
the parameter a. This fact can be used to measure the interaction J. In Figure 
3 we plot the dependence of the probability on J. The attenuation of the second 
maximum can be achieved by increasing the ratio c/a. 



0.8 
0.6 
0.4 
0.2 



5 10 15 20 

Figure 3 - Probability as a function of the interaction J for the c = 30, t — Ins and 

u = 15 GHz. 
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4 Final remarks 



We have described a method to construct some classes of exact solutions of the 
one-spin equation and the respective external field for this solution. Although 
this method can not be used to solve the one-spin equation for a specific external 
filed, it is very powerful in constructing solution whose external fields have 
some desired characteristic. As an application, we explicitly construct the exact 
solutions for some restrict in time external fields. These are a very important 
kind of fields, once the fields used in practical application usually act in a finite 
time interval. After that, we show how these results can be applied in the 
problem of two interacting spin subjected to different magnetic fields, i.e., to 
solve the two-spin equation for parallel fields. In this case, our method can be 
used to control, not only the characteristics of the external magnetic fields, but 
also the behavior of the interaction function between the spins. In a general 
manner, the exact solutions of the two-spin equation have a wide application in 
the description of two coupled quantum dots and, especially, in the construction 
of quantum gates. In order to clarify this point, we show how an arbitrary exact 
solution, not only the ones obtained with the method described here, can be used 
to obtain the operational characteristics of a Swap operation. In this analysis we 
chose the very important case of a restricted in time adiabatic variation of the 
external magnetic field, and show explicitly the dependence of the operational 
characteristics on the parameters of the external field. Besides, we describe 
the behavior of the system using explicit values of the parameters that can be 
obtained in experimental conditions. A graphical analysis of this behavior is 
presented, from where we can see that the gate can be better controlled if the 
experimental setup was assembled respecting some specific relation between the 
parameters. 
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